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Abstract 

Motivated by a recent analysis which presents explicitly the general solution, we consider the 
eigenvalue problem of the spinless Salpeter equation with a ( "hard-core amended" ) Coulomb 
interaction potential in one dimension. We prove the existence of a critical coupling constant 
(which contradicts the assertions of the previous analysis) and give analytic upper bounds on 
the energy eigenvalues. These upper bounds seem to disprove the previous explicit solution. 
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1 Introduction: The Spinless Salpeter Equation in 
One Dimension 

The spinless Salpeter equation arises either as a standard reduction of the well-known 
Bethe-Salpeter formalism [||] for the description of bound states within the framework 
of relativistic quantum field theory or as a straightforward relativistic generalization of 
the nonrelativistic Schrodinger equation. This semirelativistic equation of motion with 
a static interaction described by the Coulomb potential (originating, for instance, from 
the exchange of a massless particle between the bound-state constituents) defines what 
we call, for short, the "spinless relativistic Coulomb problem.^' 

Recently, confining the configuration space to the positive half-line (and mimicking 
thereby the effect of a "hard-core"), the relativistic Coulomb problem has been studied 
in one dimension || . This one-dimensional case may serve as a toy model which might 
prove to be instructive for the analysis of the still unsolved three-dimensional problem. 
In view of its potential importance, we re-analyze this nontrivial and delicate problem. 

The spinless Salpeter equation may be regarded as the eigenvalue equation 

H\xk) = E k \ Xk ) , k = 1,2,3,... , 
for the complete set of Hilbert-space eigenvectors \xk) an d corresponding eigenvalues 

P - (Xk\H\xk) 

k = ; i v 

\Xk\Xk) 

of a self-adjoint operator H of Hamiltonian form, consisting of a momentum- dependent 
kinetic-energy operator and a coordinate-dependent interaction-potential operator: 

H = T + V, (1) 

where T is the "square-root" operator of the relativistic kinetic energy of some particle 
of mass m and momentum p, 

T = T{p) = .J + m 2 , (2) 

and V = V(x) is an arbitrary, coordinate-dependent, static interaction potential. The 
action of the kinetic-energy operator T on an element ifi of L 2 (R), the Hilbert space of 
square- integrable functions on the real line, R, is defined by (cf. also Eq. (3) of Ref. |J) 

+oo +oo 

(Tij))(x) = — J dp J dy ^p 2 + m 2 exp[ip (x - y)] ip{y) . (3) 

— oo — oo 

In Ref. ||, the domain of H is restricted to square- integrable functions \fr(x) with 
support on the positive real line, R + , only, vanishing at x = (cf. Eq. (28) of Ref. 0): 

ijr(x) = for x < . 

This restriction may be interpreted as due to the presence of a "hard-core" interaction 
potential effective for x < 0. For x > 0, the interaction potential V is chosen to be of 
Coulomb type, its strength parametrized by a positive coupling constant a, i.e., a > 0: 

V(x) = Vc(x) = for x > . 

x 

Let the Coulomb-type semirelativistic Hamiltonian Hq be the operator defined in this 
way. 

1 The present state-of-the-art of the three-dimensional relativistic Coulomb problem 
has been reviewed, for instance, in Refs. ||, p|, ||. 
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2 Concerns — Dark Clouds Appear at the Horizon 

Now, according to the analysis of Ref . |J , the point spectrum of the Hamiltonian Hq 
consists of the set of eigenvalues (cf. Eq. (33) of Ref. ||) 

771 

E n = , 77 = 1,2,3,... . (4) 

/ a 2 

The corresponding eigenfunctions ^ n (x) must be of the form (cf. Eq. (28) of Ref. ||) 

#n(aO = j) n {x) Q(X) , 77=1,2,3,..., (5) 

where 0(x) denotes the Heaviside step function, defined here by 

Q(x) = 1 for x > , 
Q(x) = for x < . 

In particular, the (not normalized) eigenfunctions ip n (x), n = 1, 2, 3, corresponding to 
the lowest energy eigenvalues E n are explicitly given by (cf. Eqs. (37)-(40) of Ref. ||) 

if)i(x) = xexp(—/3ix) , 

( 777^ \ 
x - -g2j^j exp(-/3 2 x) , 

lf\ I 2 W ^ 3 777^1 + 7772) \ 

= x I x - -^p- x + 2 ^ ^ I exp(-/3 3 x) , (6) 

with (cf. Eq. (32) of Ref. §) 



ma a ~ 

Pn = 1 = —E n , 77=1,2,3, 

a 2 n 



n f + V 2 

and the abbreviation (cf. Eq. (26) of Ref. [|J) 



S n = ^m 2 -P 2 , 77 = 1,2,3,.... 

However, there are some facts which cause severe doubts about the validity of this 
solution: 

Boundedness from below: For coupling constants a larger than some critical value 
a c (which has yet to be determined), the operator Hq is not bounded from below. 
This may be seen, for instance, already from the expectation value of Hq with 
respect to the (normalized) trial state |$) defined by the configuration-space trial 
function 

$(x) = (p(x) Q(x) 

with 

<p(x) = 2 /i 3 / 2 x exp( — n x) , fi > , 
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and satisfying the normalization condition 

oo 

|||$)f = ($|$) = J da;|^(x)| 2 = l . 
o 

Apart from the arbitrariness of the variational parameter this trial function $ 
coincides, in fact, with the ground-state solution as given in Eqs. (|5|), (|5|). The 
expectation value of the Coulomb interaction-potential operator Vc with respect 
to the trial state |$) reads 



oo 

($|Vc|$) = -a ( dx -\ip(x)\ 2 = -a 
J x 



a 



There is a trivial (but nevertheless fundamental) inequality for the expectation 
values of a self-adjoint (but otherwise arbitrary) operator O = and its square, 
taken with respect to an arbitrary Hilbert-space state \tjj) in the domain T>(0) of 
this operator O: 



ICH '" ;! for aU G V(0) 



Application of this inequality to the kinetic-energy operator T of Eq. @) allows 
to get rid of the troublesome square-root operator: 



($|T|$) < ^/($|T 2 |$) = ^($|p 2 |$) + m 2 . 
The expectation value of p 2 required here reads 

($b 2 |$) =tf ■ 

Thus the expectation value of the Coulomb-like semirelativistic Hamiltonian Hq 
with respect to the trial state |<&) is bounded from above by 



($|#c|$> = ($\T + V C \®) < ^ii 2 + m 2 -fia . (7) 

When inspecting this inequality in the limit of large //, that is, for /jl — > oo, one 
realizes that, for a large enough, the operator He is not bounded from below. In 
fact, the expectation value of the kinetic-energy operator T with respect to the 
trial state |$), 



+ r° 4 u 3 °r \V 2 + m 2 

($|T|$)= j dx&(x)(T*)(x) = -?-j d P J P , (8) 

-oo 71 ^ " ' 

is simple enough to be investigated explicitly. For // m, this expectation value 
simplifies to 

($|T|$) = — for /j, > m . 

TV 

Consequently, in the (ultrarelativistic) limit \i — > oo, the expectation value of Hq 
behaves like 

(glgclg) 2 
hm = a . 
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This clearly indicates that for the Hamiltonian Hq to be bounded from below the 
Coulomb coupling constant a has to be bounded from above by the critical value 



2 

« c < - • 

7T 



(This upper bound on a c is, in fact, identical to the critical coupling constant a c 
found in the case of the three-dimensional spinless relativistic Coulomb problem 

00 

Upper bound on lowest eigenvalue: As rather trivial consequence of the famous 
minimum-maximum principle J7|, the expectation value 

U\H\rb) 



of a self-adjoint operator H bounded from below, with respect to some arbitrary 
state \ip) in the domain of H, D(H), is always larger than or equal to the lowest 
eigenvalue E\ of H f\ 

£,<^M for an M 6 V(H) . 

Accordingly, minimizing the expression on the right-hand side of inequality (|7|) 
with respect to the variational parameter \i yields a simple analytic upper bound 
Ei on the ground-state energy eigenvalue E\ of the Coulomb-like semirelativistic 
Hamiltonian He' 

Ei < Ei 

with 

Ei = m Vl - a 2 . (9) 

The same analytic upper bound on the ground-state energy Ei has been found in 
the case of the three-dimensional spinless relativistic Coulomb problem 
Reality of this latter expression requires again the existence of a critical coupling 
constant a c and indicates that this critical value of a is less than or equal to 1: 

a c < 1 . 

Moreover, at least for the energy eigenvalue Ei corresponding to the ground state 
of the Hamiltonian He, the supposedly exact value of Eq. (Q), 

Trb < (10) 

is in clear conflict with the naive upper bound Ei of Eq. (|]): 

Ei 
Ei 

and therefore 



a 



l 



Ei < Ei for a > 



This statement constitutes what is sometimes simply called "Rayleigh's principle." 
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For larger values of the Coulomb coupling constant a, the upper bound (§) on the 
ground-state energy can be easily improved by fixing in the expectation value @ 
of the kinetic-energy operator T the variational parameter \x to the value fj, = m. 
In this case, this expectation value reads 

($|T|$) '"' 



7T 



Accordingly, the ground-state energy eigenvalue E 1 is bounded from above by 

"4 



Ei < y n | /.>/ 

For the Coulomb coupling constant a in the range 

2 
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7T 



14 2 
< a < - , 



7T 
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the above expression represents a genuine improvement of the upper bound (f|). 

Eigenstate expectation values vs. eigenvalues: The expectation value @ of the 
kinetic-energy operator T with respect to the trial state |<&) may be written down 
explicitly: 

o arccos 
($|T|$) = -m + 



7T 



m 



1 



ill 

m 2 



Now, for 



ma 



the trial function $ coincides with the normalized ground-state eigenfunction y&i. 
In this case, the corresponding expectation value of the Hamiltonian Hq becomes 

2 

— (a + (1 + a 2 ) arccot a) — a 2 



(*i|#c|*i) 



rn 



a- 
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(12) 



Unfortunately, the above expectation value does not agree with the ground-state 
energy ([TOf) deduced from Eq. (|): 

Orthogonality of eigenstates: Eigenstates \Xi), i = 1, 2, 3, . . . , of some self-adjoint 
operator H corresponding to distinct eigenvalues of H are mutually orthogonal: 

(Xi\Xk) oc 8* , h k = 1,2,3, . . . . 
This feature is definitely not exhibited by the overlaps 



dx*i(x)V k (x)= / dxtf(x)ip k (x) , z,fc = l,2,3, 



(*i|* fc > 



of the lowest eigenfunctions ^(x), z = 1, 2, 3, given in Eqs. (|6|). For instance, 
the overlap (^il^) of the ground state and the first excitation j^) is given 

by 

2[3£ 2 &-m 2 (A+&)] 
revealing thus, beyond doubt, the non-orthogonality of the vectors and l^)- 
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3 Exact Analytic Upper Bounds on Energy Levels 

In view of the above, let us try to collect unambiguous results for the one- dimensional 
spinless relativistic Coulomb problem. With the help of the definition (|3|) of the action 
of a momentum-dependent operator in coordinate space, it is easy to convince oneself 
of the validity of the operator inequality 

p 2 

T < T NR = m + 



2m ' 

the relativistic kinetic-energy operator T is bounded from above by its nonrelativistic 
counterpart Tnr: when introducing the Fourier transform ip(p) of the coordinate-space 
representation ip(x) of the Hilbert-space vector 

+oo 

Tp(p) = — / dx exp[— ipx] ip(x) , 



'2tx 



-oo 



one finds 

+oo 

(V|T NR -r|v> = / dif(i)[W)(i)-(^)(i)] 



-oo 



dp \ifj(p)\ 2 ( m + — J p 2 + m? ] 

V 2m v / 

— oo v ' 

> o . 

Hence, adding the Coulomb interaction potential Vc, the semirelativistic Hamiltonian 
He is, of course, bounded from above by the corresponding nonrelativistic Hamiltonian 

-^c,nr: 

He < #C,NR = 2nr + Vc ■ 

Now, upon invoking the minimum-maximum principle (which requires the operator 
Hq to be both self-adjoint and bounded from below) and combining this principle with 
the above operator inequality, we infer that every eigenvalue E n , n — 1, 2, 3, . . . , of H c 
is bounded from above by a corresponding eigenvalue E n>NR , n = l,2,3,...,of i?c,NR : Q 

E n < -E„,nr for n = 1, 2, 3, . . . . 
It is a simple and straightforward exercise to calculate the latter set of eigenvalues: 

( a2 \ 
#n,NR = m 1 1 - — -J I , n = 1, 2, 3, . . . . 

These upper bounds on the energy eigenvalues E n may be easily improved by the same 
reasoning as before. Introducing an arbitrary real parameter r\ (with the dimension of 
mass), we find a set of operator inequalities for the kinetic energy T [§, f|, [10| , namely, 



T < p2 + m2 + r l 2 forall,>0 
2?y ' 



3 The line of arguments leading to the general form of this statement may be found, 
for instance, in Refs. || |J. It is summarized in Appendix [A]. For a rather brief account 
of the application of these ideas to the three-dimensional spinless relativistic Coulomb 
problem, see, e.g., Ref. [§]. 
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and, consequently, a set of operator inequalities for the Coulomb-type semirelativistic 
Hamiltonian H c H ||: 

^ _|_ jy)^ _|_ ^2 

#c < #cfa) = + Vb for all r] > . 

2 Tj 

Accordingly, every eigenvalue E n , n = 1, 2, 3, . . . , of H c is bounded from above by the 
minimum, with respect to the mass parameter rj, of the corresponding eigenvalue [§, [| 

„2 S 



2 



m 2 + ?7 2 I 1 - 



q- 



n = 1,2,3, 



of H c {rf): 



a 2 



E n < min E n c(ri) — m \ 1 for all a < a c 

?7>o ' y 



For n = 1, this (variational) upper bound coincides with the previous upper bound (|9|). 
It goes without saying that these upper bounds are violated by the energy eigenvalues 
E n given in Eq. M): 



means 



1 a 4 

min E n c (r]) — \ 1 7 < 1 for a 7^ , for all n = 1, 2, 3, . . . 

^>o ' V n 4 

min .E n c (r]) < E n for a 7^ , for all n — 1, 2, 3, ... ! 

r?>0 ' 



Moreover, for \x = ma, our generic trial state |$) becomes the lowest eigenstate of 
the nonrelativistic Hamiltonian if CNR , corresponding to the ground-state eigenvalue^ 

( « 2 \ 
E xm = m I 1 - — 1 , 

which may be easily seen: 

(T NR v?)(x) = ( m - -^-r J p(ar) = ( m - ^- + — - ) <p(x) for a; > 
\ 2max z J \ 2 m mxj 

implies (with fi = ma) 

Hc,nr\&) = -Ei.nrI^?) . 

It appears rather unlikely that the same functional form represents also the eigenstate 
of the semirelativistic Hamiltonian He- 



4 Summary, Further Considerations, Conclusions 

This work is devoted to the study of the one- dimensional spinless relativistic Coulomb 
problem on the positive half-line. Assuming a (dense) domain in L 2 (R + ) such that the 
semirelativistic Coulombic Hamiltonian Hq defined in the Introduction is self-adjoint, 
analytic upper bounds on the energy eigenvalues E k , k = 1,2,3, ... , have been derived: 



E k <m]ll-— for all k= 1,2,3, ... . (13) 

Surprisingly, the explicit solution presented in Ref. || does not fit into these bounds. 

4 The Coulomb problem involves no dimensional parameter other than the particle 
mass m. Therefore, both the energy eigenvalues E n and the parameter(s) [i have to be 
proportional to m. 
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In order to cast some light into this confusing situation, let us inspect the action (|3|) 
of the kinetic-energy operator T in more detail. Consider not normalized Hilbert-space 
vectors \<& n ), n = 0, 1, 2, . . . , defined, as usual, by the coordinate-space representation 

$ n (x) = x n exp(-fix) Q(x) , fi > , n — 0,1,2,.... 

These vectors certainly belong to the Hilbert space L 2 (R) for all n = 0, 1, 2, ... , since 

|||$n)|| 2 = ($n|$n> = / &X \<$> n (x)\ 2 = f dxX 2 " exp(~2 fJL x) = [ ™ / < OO . 

-oo W 

However: The norm ||T|$„)|| of the vectors T|<E> n ), n — 0, 1, 2, . . . , may be found from 

||T|*„>||*= Td,| (r *„)wr = M^Td P . > f + m ' 



2n J 1 (p 2 + n*)** 1 ' 

— ^ — oo 

This observation might be a hint that the vector |$o), that is, &o(x) = exp(— \x x) G(x), 
does not belong to the domain of the kinetic-energy operator T. If this is indeed true, 
it is by no means obvious how to make sense of Eq. (16) of Ref. || for the case n = 0. 

Trivially, if Eq. (16) of Ref. || is correct for n = 0, all these relations for arbitrary 
n = 1,2,... may be obtained by a simple differentiation of the relation for n = with 
respect to the (generic) parameter fi, taking advantage of 




T x n exp(— (jlx) = . 

Similarly, it is somewhat hard to believe that Eq. (16) of Ref. [f| holds for n = 1. In 
our notation, Eq. (16) of Ref. would read for n = 1 



Crooo*) = (s + ^) 



with 

S = \Jm 2 — n 2 . 

Considering merely the norms of the vectors on both sides of this equation, we find, for 
the norm of the vector on the left-hand side, 

im»oii' = =^ 

but, for the norm of the vector on the right-hand side, 

9 4 i 4 



m + fj, 



Afi 3 S 2 ' 

These two expressions for the norms become equal only for the — excluded — case \i = 0. 
Unfortunately, precisely the above relation forms the basis for the assertion in Ref. 
that $i (x) with [i = (5\ is the ground-state eigenfunction of the ( "hard-core amended" ) 
one-dimensional spinless relativistic Coulomb problem as defined in the Introduction. 

In conclusion, let us summarize our point of view as follows: The energy eigenvalues 
Ef., k — 1, 2, 3, . . . , of the one-dimensional spinless relativistic Coulomb problem (with 
hard-core interaction on the nonpositive real line) are bounded from above by Eq. ( JT3| ) . 
For the ground-state energy eigenvalue Ex, this upper bound may be improved to some 
extent, by considering appropriately the minimum of the bounds of Eq. (|TT|), Eq. ([12]), 
or Eq. ( |i"3"|) for k — 1, that is, Eq. (||). To our knowledge, these upper bounds represent 
the only information available at present about the exact location of the energy levels 
of the ( "hard-core amended" ) one-dimensional spinless relativistic Coulomb problem. 
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A Combining Minimum— Maximum Principle with 
Operator Inequalities [9], gj] 

There exist several equivalent formulations of the well-known "min-max principle" 0] . 
For practical purposes, the most convenient one is perhaps the following: 

• Let H be a self-adjoint operator bounded from below. 

• Let Ek, k = 1,2,3, ... , denote the eigenvalues of H, defined by 

H\Xk) = E k \xk) , k = 1,2,3, ... , 
and ordered according to 

E\ < E2 ^ -E3 ^ • • • • 

• Consider only the eigenvalues E k below the onset of the essential spectrum of H . 

• Let Dd be some d- dimensional subspace of the domain T>(H) of H: C T>(H). 
Then the fcth eigenvalue Ek (when counting multiplicity) of H satisfies the inequality 

E t < 8 up for fc= 1,2,3,.... 

\4>)€D k (W> 

The min-max principle may be employed in order to compare eigenvalues of operators: 

• Assume the validity of a generic operator inequality of the form 

H < O . 

Then 

P _ (xk\H\ Xk ) , (VW) . W>W> 

E k = — : — < SUp — 7777T— < SUp 



• Assume that the A;-dimensional subspace Dk in this inequality is spanned by the 
first k eigenvectors of the operator O, that is, by precisely those eigenvectors of O 
that correspond to the first k eigenvalues E±, E 2 , . . . , Ek of O if the eigenvalues of 
O are ordered according to 

E\ ^ E2 ^ i?3 ^ ■ ■ ■ • 

Then 

sup = E k . 

WeD k {ip\ip} 

Consequently, every eigenvalue E k , k = 1, 2, 3, . . . , of H is bounded from above by the 
corresponding eigenvalue E k , k = 1, 2, 3, . . . , of O: 

E k < E k for k = 1,2,3,... . 
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